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We have searched for the Lepton Flavor Violating decay r~ fi'rj using a data sample of 84.3 
fb~^ accumulated with the Belle detector at KEK. The ry-meson was detected through the decay 
modes: — > 77 and 7r+7r~7r°. No signal candidates are found, and we obtain an upper limit for the 
branching fraction B{t^ l^ v) < 3.4 x 10"^ at the 90% confidence level. 

PACS numbers: 11.30.-j, 12.60.-i, 13.35.Dx, 14.60.Fg 



Among the possible Lepton Flavor Violating (LFV) 
decays of the r-lepton, r~ — s- ^~ri is the process that 
provides the most stringent bound on Higgs-mediated 
LFV. Sher has pointed out that a flavor non-diagonal 
lepton-lepton-Higgs Yukawa coupling could be induced if 
slepton mixing is large. The ^-r-Higgs vertex is particu- 
larly promising since mixing between left-handed smuons 
and staus is large in many supersymmetric models 0. 
This mechanism initially led various authors to study 
the enhancement of the LFV decay t — > 3/i in the mini- 
mal supersymmetric standard model (MSSM). However, 
Sher's results indicate that fi^rj is enhanced by a 

factor of 8.4 compared to r — ^ 3/^, due mainly to a color 
factor and the mass-squared dependent Higgs coupling 
at the Higss-s-sbar vertex. In some models with rea- 
sonable assumptions about MSSM parameters [1,3] the 
T~ — > n~r] branching fraction is given by 

N ^o. 6 /tan/3\VlOO GeVV 

where itia is the pseudoscalar Higgs mass and tan/3 is 
the ratio of the vacuum expectation values {{Hy) / {Hd))- 
In such models, — > and t 3fj, are particularly 
sensitive to LFV at large tan/3. 

Previous experimental studies of r~ ii~rj by AR- 
GUS and CLEO ^ set 90% confidence level upper 
limits on the branching fraction of 7.3 x 10~^ from 0.387 
fb-i of data, and 9.6 x 10"^ from 4.68 fb"! of data, re- 
spectively. We present here a new search based on a data 
sample of 84.3 fb~^, equivalent to 76. 9M t+t^ pairs, col- 
lected at the T{AS) resonance with the Belle detector at 



the KEKB asymmetric e+e" colliderjg. A description 
of the detector can be found in Ref. [7|. 

For Monte Carlo (MC) studies, the following pro- 
grams have been used to generate background (BG) 
events: KORALB/TAUOL A for t+t" processes, 
QQ 13 for BB and continuum, BHLUMI [13 for Bhabha, 
KKMC 113 for MM and AAFH [H for two-photon pro- 
cesses. The T~ fi~r] decay is initially assumed to have 
a uniform angular distribution in the r's rest frame. The 
Belle detector response is simulated by a GEANT3 
based program. Most kinematical variables are evalu- 
ated in the laboratory frame, unless denoted by the su- 
perscript "CM" in which case they are evaluated in the 
center-of-mass frame. Two 77 decay modes are consid- 
ered in this analysis: 77 ^ 77 (6 = 39.43 ± 0.26%) and 
r] TT+TT^Tr" {B = 22.6 ± 0.4%) 

For rj ^ ^j, we search for events containing exactly 
two oppositely charged tracks and two or more photons, 
two of which form an rj. The events should be consistent 
with a T+T" event, in which one r decays to iirj and the 
other T decays to a charged particle other than a muon 
with any number of 7's and neutrinos. 

To select candidate events we require the momentum 
of each track, p, and the energy of each photon, E^, to 
satisfy p > 0.1 GcV/c and > 0.1 GcV. The tracks 
and photons are required to be detected in the barrel or 
endcap of Belle: -0.866 < cos6' < 0.956. To exclude 
Bhabha, fifj, and two-photon events, the total energy is 
constrained between 5 and 10 GeV in the CM frame, as 
shown in Fig. ^a). 

In the CM frame the events are subdivided into two 
hemispheres by a plane perpendicular to the thrust axis. 
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FIG. 1: Some kinematical distributions from MC and data: a) 
total energy, -E^taii b) cosine of the opening angle between the 
fj, and 77 on the signal side, c) m^igg vs. pmiss for signal MC 
events, d) m^i^g vs. Pmiss for generic decays of t'^t' MC. The 
shaded histograms represent the signal r"'"r~ MC, the hatched 
histograms represent the combined i-i'^ fJ.^ and Bhabha MC, 
the open histograms represent the combined BB, continuum 
and two-photon MC and the closed circles represent the data. 
The selected regions in a), b) are indicated by the arrows. 
The selected region in c) , d) is the area between the two lines: 



Pv 

0.869. 



> -2.615 X mi, 



0.191 and 



> 1.238 X mi, 



The signal side should contain a muon and two pho- 
tons. The muon is identified as a track having a fi 
probability > 0.9 |l5| . An rj meson produced in 
a two-body t decay has on average a higher momen- 
tum than rj mesons from other sources. Therefore, a 
photon from rj decay is required to have a rather high 
energy Ej > 0.22 GeV. To reduce background, events 
are rejected when two 7's, one from the signal side and 
the other from the tagging side (7'), have a resolution- 
normalized Tr'^-mass in the range —5 < SJ^^, < 5, where 
Sf^, = (m^y - 0.135 GeV/c2)/a;;^, and a;!^, is in the 
range 5-8 MeV/c^. This Tr° veto rejects 86% of BG events 
while retaining 75% of the signal events. To further re- 
duce BG events, the cosine of the opening angle between 
the ^ and 77 on the signal side is required to satisfy 
0.5 < cos0p_^-y < 0.95, as shown in Fig. H^b). 

In the tagging side hemisphere, the charged track 
should not be a muon {V^ < 0.6 is imposed), but may be 
either an electron, i.e. have an e probability Ve > 0.9 0| 
or a hadron {V^ < 0.6 and Ve < 0.9). If an electron is 
found, the number of photons and the electron momen- 
tum are constrained by n.^ < 2 and pe > 0.7 GcV/c. If a 
hadron is found, the constraints are > and phad > 
0.1 GeV/c. 



The following two criteria are imposed on the miss- 
ing momentum and energy in the event. To ensure that 
the missing particles are neutrinos rather than 7's or 
charged particles that fall outside of the detector accep- 
tance, we require that the direction of the missing mo- 
mentum should satisfy —0.866 < cos^miss < 0.956. Be- 
cause neutrinos are emitted only on the tagging side, the 
direction of the missing momentum should be contained 
on the tagging side: cos ^thrust- miss < -0.55. The cor- 
relation between the missing momentum, Pmiss, and the 
missing mass squared, TO^jgg, shown in Figs.^c) and (d) 
for signal and generic t~^t~ MC, is utilized for additional 
BG rejection. 

The 77 candidate is selected based on the signal-side 
77 invariant mass in terms of the resolution-normalized 



77-mass, 



-5 < 



77 



< 3, where S'lJ^ 



(™7 



0. 547 GeV/c2)/(7^^ and a!^^ is 12 MeV/c^. The resuhing 
Sl^^ distributions for signal and generic r+r^ MC and 
data are shown in Fig. |2Ia). 

The application of these selection criteria to the data 
set results in a total yield of 18 events. The detection 
efhciency is measured from MC studies to be £(27) = 
9.3%. In MC, small backgrounds from the three follow- 
ing processes survive: 8.6±2.2 events from generic r+r~, 
2.5±1.8 events from /i/i and 5.8±2.2 events from the con- 
tinuum. 

For the rj — > Tr+Tr^Tr*^ mode, we search for events con- 
taining four charged tracks (net charge = 0) and two 
or more photons. Because of the higher multiplicity 
compared to the 77 —^ 77 mode the detection efficiency 
is smaller; however, the extra reconstruction constraint 
in the 77 decay chain improves the background rejection 
power. The selection criteria are similar to those in the 
rj — > 77 case with the differences listed below. 

The minimum photon energy is reduced from 0.1 GeV 
to 0.05 GeV, since the photons from this decay mode have 
a softer energy distribution compared to those in 77 — > 
77. The signal side hemisphere should have three tracks 
and two or more photons. One track must be a muon 
{Vfi > 0.9), but particle identification is not performed 
on the other two tracks — they are treated as pions. 
We also require that one 7r° be reconstructed from the 
photons in the signal hemisphere, such that — 5 < S^^ < 
5. Figure[2Ib) shows the reconstructed mass of 77. 

After the cuts, 67 events remain in the data, while 
the generic r+r" MC predicts a contribution of 38.0±4.6 
events, and the continuum MC predicts 15.6 ± 3.5 events. 
The detection efficiency is e(37r) = 5.6%. 

The final evaluation of the number of signal candidates 
is performed by defining a signal-region in the M^^-Ai? 
plane, where the candidate /ir] system should have an 
invariant mass (M^ri) close to the r-lepton mass and 
an energy close to the beam-energy in the CM frame, 

1. e. AE ^ £;CM _ £;™^ c± 0. Figures E^a) and 
IHl^b) show scatterplots of the signal MC in the M^^)- 
AE plane for the 77 — > 77 and rj — > 7r"'"7r~7r'' modes. 
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FIG. 2: (a) invariant mass of 77 in terms of the resolution 
normalized 7;-mass, S^^, in the 77 ^ 77 case, and (b) 77-mass 
from 77 — > 7r^7r"7r'' reconstruction. Signal and generic r^r^ 
MC distributions are indicated by the shaded and open his- 
tograms, respectively. The selection region is indicated by the 
arrows. 



respectively. The signal exhibits a long low-energy tail 
due to initial-state radiation and calorimeter energy leak- 
age for photons. By reproducing the M^^ and A£' dis- 
tributions around the peak with an asymmetric Gaus- 
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sian function, the il/^^ and AE resolutions are evalu- FIG. 3: Final event distributions over a ±10cr region in the 



ated to be a 



low/ high 



= 25.8±0.7/15.3±0.4 MeV/c^ and 



low/high 



69.7±3.0/34.7±1.2 MeV for the 77 



77 

mode, and cr^°J'/^'sb^ 13.8±0.4/9.0±0.4 MeV/c^ and 



low/high _ 



44.4±2.3/22.6±1.3 MeV for the 77 tt^ 



mode, where the "low/high" superscript indicates the 
lower/higher energy side of the peak. To optimize the 
sensitivity, we take an elliptically shaped signal region in 
the M^rj-AE plane, with a signal acceptance of 17=90%, 
as shown in Fig. 

Figure 13 shows the final data distributions for a 
±10(Tjvfj,^/A_E region in the M^j^-Ai? plane. In the signal 
region, there are no events in either the data or back- 
ground MC. Outside the signal region, 7 events for the 
77 — > 77 mode and 2 events for the 77 tt+tt^tt'^ mode 
are observed in data, while MC predicts 3.7 ± 2.4 and 
O.Oto^o events, respectively. The observed data yields arc 
consistent with MC. The BG yield in the signal region, 
estimated from the sidebands, is found to be 0.5 ±0.2 for 
77 — > 77 and 0.0^0 events for 77 7r"''7r~7r°. 

As no events are observed, an upper limit on the num- 
ber of events is set using a Bayesian approach, which 
gives So = 2.3 at 90% C.L. The upper limit on the branch- 
ing fraction, at 90% C.L., is given by 



M^r)-AE plane for (a) the rj —> ^"f and (b) 77 — > n TT~n 
modes. The ellipses are the signal regions with an acceptance 
of n = 90%. The data are indicated by the open circles, and 
the signal MC events are plotted as dots. 



The systematic uncertainties on the detection sensitiv- 
ity, 2{e^l X Brf) X N^+T- , arise from the track reconstruc- 
tion efficiency (2.0% in the 77 -f-f mode and 2.0% in the 
77 — > TT+TT^TT^ mode), 77 reconstruction efficiency (2.0% 
and 4.2%, the latter of which includes the uncertainties of 



tracking efficiency for 77 



^ttO), TfO veto (5.5% and 



B{t v)< 



So 



2 {efl X B^) X 7V^4 



(2) 



where S,, is the branching fraction of 77-decay to either 
77 or 7r+7r^7r°. The calculated upper limits, at 90% C.L., 
are thus found to be 4.6 x 10~^ for the 77 — > 77 mode, 
and 13.1 x lO^*" for the 77 tt+tt^tt*' mode. Combining 
the two decay modes, we obtain eil x S,, ~ 4.4% and 
B{t- fi-ri) < 3.4 X 10-^ at 90% C.L. 



none), muon identification efficiency (4.0% and 4.0%), 
trigger efficiency (1.4% and 1.4%), beam background 
(2.3% and 2.1%), luminosity (1.4% and 1.4%), 6^ (0.7% 
and 1.8%) and MC statistics (1.3% and 2.1%). Adding 
all of these components in quadrature, the total uncer- 
tainty is evaluated to be 8.1% for 77 — > 77 and 7.3% for 
77 — > 7r"'"7r~7r''. For the combination of the two decay 
modes the systematic uncertainty is ±7.9%. 

This systematic uncertainty is included in the upper 
limit following Ref. [l^ . where the detection sensitiv- 
ity, 2(eil X BtiJNtt, is modelled by a Gaussian distribu- 
tion having a width given by the systematic error quoted 
above. There is no appreciable effect on the branching 
fraction, B. 

The angular distribution of the ^^rj decay has a 

strong dependence on the LFV interaction structure 
and spin correlations between the r's at the signal and 
tagged sides must be considered. To evaluate the max- 
imum possible variation, V— A and V+A interactions 
are assumed; no statistically significant difference in the 
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FIG. 4: Experimentally excluded ruA —tan l3 parameter space. 
The result of this experiment using Q is indicated by the 
shaded region together with the regions excluded by LEP ^ 
and the Tevatron ITllll. 



Mfj,rj-A.E distribution or in the efficiency is found com- 
pared to the case of the uniform distribution. 

As a result, we obtain an upper limit on the branching 
fraction for the Lepton Flavor Violating — > fi^rj decay 
of 



fi^r]) < 3.4 X 10" 



(3) 



at 90% C.L. This result improves the previous upper 
limit, B{t- fi-j]) < 9.6 X 10"^ 0, by a factor of 
30. 

Using Eq.(^, which was derived in a seesaw MSSM 
with a specific neutrino mass texture, our upper limit 
restricts the allowed parameter space for ttia and tan/3, 
as indicated in Fig. ^ where our boundary is indicated 
in the cases of 90% and 95% C.L. Figure^ also shows 
the 95% C.L. constraints from high energy collider ex- 
periments at LEP and CDF [l^. Our result has a 
sensitivity close to that of the CDF experiment, achieved 
by searching for pp A/(j)bb — > bbbb, where </> is a CP- 
even neutral Higgs state and A is a CP-odd state in the 
MSSM. 
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